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STABILITY OF BEAM COLUMNS ABOVE 
THE ELASTIC LIMIT 


Robert L. Ketter,” J.M., ASCE 


SYNOPSIS 


This paper, which discusses the behavior of pin-ended, steel, wide-flange 
beam columns, is the development of a virtual displacement method of solu- 
tion to the inelastic instability problem. The method is equally applicable to 
the solution of elastic or inelastic beam-column problems but its greatest ap- 
plication seems to be in the inelastic field. Using the proposed method of 
solution, the critical value of loading for members subjected to concentric, 
eccentric and or lateral loads, can be determined in a fraction of the time 
required by existing solutions. 


INTRODUCTION 


The problem under consideration is the determination of the maximum 
loading a beam-column type member can carry. The importance of the prob- 
lem is evident, since if it is possible to predict the "collapse" loading of a 
member, it is also possible to evaluate for a member with a given loading the 
true factor of safety. 

A large number of variables enter in the solution of problems of this type, 
and to reduce these in the development that will follow, it will be assumed 
that: 

a. The mode of failure is that of bending in the plane of the 

applied moments. Furthermore, this plane is coincident with a 

principal axis of the section, (Lateral-torsional and/or local buckling 

are not considered.) 

b, The material under study possesses stress-strain properties as 

shown in Fig. 1. 

c., Members are originally straight, free from accidental eccen- 
tricities, and of uniform cross-section along their lengths. 


To illustrate the phenomena being studied, consider the typical load- 
deflection curves of Figs. (2) and (3). 

For the axially loaded member, Fig. 2 (elastic behavior), axial thrust, P, 
can be increased from zero to a certain critical value, Pe, with the member 
remaining in the straight equilibrium position. However, on application of 
the slightest additional increment of thrust above the critical value, lateral 
deformation occurs, Thus at the instant of bending there exist two adjacent 
positions of equilibrium of the same load; that is to say, the member is in- 
different with regard to either a straight or a deflected configuration. This 
condition, known as bifurcation of the equilibrium position, is called a 
"condition of instability.” 


*Research Instructor, Fritz Eng. Lab., Lehigh Univ., Bethlehem, Pa, 
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For the eccentrically loaded column problem of Fig. 3 a different type of 
behavior is observed, First of all, the member starts to deflect on the 
slightest application of thrust, Elastically this deformation increases with 
the maximum strength approaching the Euler load, However, as deflections 
increase with load, the outermost compressive fiber stress at the center of 
the column increases as the product of P and yp. Because of this rapidly in- 
creasing stress situation, there is reached a load at which yielding of the 
member takes place, (e.g., point "a" in Fig. 3). As the load is further in- 
creased, yielding progresses across and along the member, thereby reducing 
its resistance to further loading. Finally, there is reached a point at which 
an increase in load is impossible, the member progressively deflecting as 
load remains constant. At this instant, the member is indifferent with regard 
to two possible adjacent equilibrium configurations and it becomes unstable, 

Herein then, "instability" is considered as that phenomenon where under 
a given constant load there exist two possible adjacent equilibrium positions. 

To determine this condition of instability, several basic methods of solu- 
tion are available. Three of these are as follows: 

a. If the force required to displace the loaded column into a dif- 
|, ferent, deflected position equals zero, then the member is said to be 
indifferent with regard to which of the positions it will assume, and a 
condition of instability exists, 
“< b. If the total work associated with a loaded column remains the 
same as the member is deflected, that is, the rate of change with re- 
spect to deflection equals zero, indifference again has taken place. ; 
c. If it is possible to determine the axial load - lateral deflection 
relation, the member will become unstable when 


since for that value of axial load there exist two adjacent equilibrium 
positions, 


EXISTING SOLUTIONS 


In the main, previous solutions of the inelastic instability, beam-column 
problem have been based on the criterion of instability that involves the de- 
termination of the maximum point of a load-deflection curve. In general each 
of those solutions has been based on the work of Karman, who, starting from 
an assumed strain distribution pattern, (and thus a stress distribution pat- 
tern), calculated the corresponding moment, thrust and curvature values 
which would have produced such a stress configuration. Since the three 
equations defining these variables, M, P and @, contain five unknowns (see 
Fig. 4), one expression of the form 


Mint {(¢, P) (1) 


can be obtained by eliminating the two cross-section variables, (@ and Oy), 
which define the generalized yield pattern. 

Consider now a typical beam-column (Fig. 5), Since equilibrium must be 
satisfied at all times, i.e. 


Mext = Mint (2) 
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and since the equation for external moment, Me xt at any point along the 
column can be expressed as a function of P, y, lateral loads, etc., i.e. 


Mext = & (P, y, lateral loads, etc.), 


solution to equation (1) will be of the form 
P = h (¢, y, lateral loads, etc.) (3) 


At this point the various methods of solutions differ. Karman and 
Chwalla(3) used numerical or graphical integration to solve equation (3) for 
direct values of P versus y. Jezek, (7) Westergaard and Osgood, (16) and 
others assumed a known deflection curve such that ¢ could be obtained di- 
rectly in terms of y by differentiation. Both of these methods of solution 
have as their objective the determination of a means of describing thrust as 
a function of deflection, 


P = f(y), (4) 


for which, according to criterion (c), instability occurs when 


(See Fig. 3.) 

The difficulties involved in either of these methods of solution arise from 
the fact that for any conventional type of cross-section other than a simple 
rectangle or circle, it is impossible to write a continuous function for Eq. 
(1); therefore curves or graphs must first be devised to supply the M-P- ¢ 
information corresponding to various stress conditions.(9) Proceeding from 
these plots it is necessary to integrate suitable ¢ values and obtain the actual 
load-deflection curve. 

To circumvent this difficulty, Jezek(7) introduced in his approximate solu- 
tion a "shape factor"* the purpose of which was to relate the relative stiff- 
ness of the section in question to that of the rectangle. This type of solution 
lends itself readily to the formulation of design equations* which are com- 
paratively easy to solve. 

One notably different type of solution to the eccentric instability problem 
was presented by Ros, (11), (12) who reasoned that if the rate of increase in 
effective, internal, resistance eccentricity, M/P, is the same as the rate of 
increase in the external applied centerline eccentricity, then instability will 
occur. Using this condition he obtained a graphical solution for the rec- 
tangular section problem by expressing the internal eccentricity as a function 
of the sum of the outer fiber strains (a measure of curvature). 

It has been shown that the presence of residual stress is of definite in- 
fluence on the strength of columns. (6), (9) This makes solution by the above 
methods even more complicated. 

To overcome the difficulties encountered in each of these solutions, the 
method presented in the next section of this paper was developed. The me- 
chanics are similar to those presented by Ros 12) for the rectangular 
section. 


*Page 45, Ref. (1). 
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THE PROPOSED VIRTUAL DISPLACEMENT 
METHOD OF SOLUTION 


Bending instability can be considered as a breakdown in the resistance the 
member offers to further bending. Since bending is synonymous with mo- 
ment application, it should be possible, then, to determine for a given loading 
condition if the section offers sufficient resistance to further bending by dis- 
placing the loaded specimen a virtual amount and noting the change in exter - 
nal applied moment at the most critical section and the internal moment of 
resistance offered at that same section. If the increase in external moment 
is greater than that which the internal stress pattern can supply, the member 
fails since resisting moment cannot "keep up" with applied moments. When 
these two ratios of increase are the same a condition of indifference exists. 
The necessary condition for instability can then be stated: 

When external moment at the most critical section increases at a rate 

equal to the rate of increase in resistance moment supplied by the in- 

ternal stress system at that same section, the member becomes un- 
stable. 
Writing in the form of an equation 


4M 
Ag 


where AM/A dis the rate of change in moment with respect to curvature,¢. 

For illustration consider the axially loaded column shown in Fig. 2 as- 
sumed, for ease of solution, to buckle in the elastic stress range. As the 
thrust, P, is increased the potential of the external system to subject the 
various cross-sections along the member to bending moment in addition to 
applied thrust increase as the thrust, i.e. 


AMex; = Pay (6) 


. 4M 


Ad (5) 


ext int 


Whereas, internally (as shown in Fig. 6) the resistance to this possible ex- 
ternal moment is constant at a value defined by the equation 


AMint = ELA (7) 
AMext and AMjnt represent changes in moment when going from the straight 
(original) to the deflected position (a virtual displacement away). 
Assuming the deflection curve is of the form 


7x 
Ay = Ay, cos i? 


then 
62 n2 
4g, = —> cos —) = Ay, 
or 
Ay. =4¢ L? (8) 
Yo Oo 


where the zero subscript denotes conditions at the center of the column (the 
most critically deformed section). Substituting Eq. 8 in Eq. 6, 
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Writing this expression and Eq. (7) as the change in moment divided by the 
virtual curvature to obtain the rate of change in moment with respect to Do 
gives 


Each of Eqs. 9 can be plotted as function of P and (AM/A@), (see Fig. 7). 
There is reached a point at which the rate of external potentiality of moment 
application (4M/A¢ |e ) exactly equals the rate of internal resistance to lat- 
eral deformation (AM/A@ |;,4)- At that point the member can assume either 
the straight or the deflected position, and therefore a condition of instability 
exists. Above that critical value of load, P,,, it is impossible for the inter - 
nal change in resistance moment to "keep up” with the potential change in ex- 
ternally applied moments and therefore the member collapses. To evaluate 
the magnitude of the load corresponding to this critical condition, use the in- 
stability condition of Eq. 5, 


4M,; 4M, 
a 


which when solved for P,, gives 


Por 4 L2 
the Euler buckling load for a straight, axially loaded, elastic column. 

Had the thrust been applied eccentrically and/or a lateral load super- 
imposed thereby introducing bending from the start, the following would have 
resulted, In the original deflected equilibrium position (see Fig. 8), prior to 
subjecting the member to the virtual disturbance, 


M P( ) RL 
= P(e + yo) + — 
Olext 4 


whereas, after the Ay, change 


L 
M, gt 


4M, L2 
= P — 
(9) 
0 
| 
or 
2 
L 
EI, 
cr 2 
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The increase in external moment caused by this disturbance is equal to the 
difference between the final and original moment values, i.e. 


= M’ 
ext ext 


ext 


4M 


s P(Ay,) 


ext 


It should be noted that this is the same as for the axially loaded column prob- 


lem shown previously.* 
Internally, again considering elastic conditions for simplification, (see 


Fig. 9), resistance to bending is constant at a value 


4M 
4¢ 


= EI 


int 


Solution to this problem then is the same as for the pure axial load case, in- 
stability occurring at the Euler buckling load. However from the discussion 
of Fig. 3 and from the bending resistance plot of Fig. 10 it is obvious that this 
condition can never be reached and satisfy the assumption of an elastic inter- 
nal stress situation. To solve the problem then, it is first necessary to 
determine internal bending stiffness above the elastic limit. 


Bending Stiffness Above the Elastic Limit 

As was discussed in Section III (Fig. 4) the moment-curvature relation 
above the elastic limit is a function of axial thrust as well as cross-sectional 
dimensions and mechanical properties of the material. Thus to be able to 
use the resulting M-¢ curves for integration it is necessary to first deter- 
mine compatible moment and thrust quantities corresponding to given % 
values. This can become a major problem for anything other than simple 
shapes; however, by modifying the problem under study to one where axial 
thrust is held constant and moments increased to collapse, a relatively sim- 
ple, semigraphical solution is available(9) whereby compatible M, P,¢ values 
can be determined for any type of cross-section having one axis of symmetry. 

Two sets of curves resulting from that previous study are here shown as 
Fig. 11 and Fig. 12 and are the M- @plots (at constant P-values) for the 
8WF31 (strong axis bending) and for the rectangular section respectively. 
These plots have been made non-dimensional by dividing the moment value by 
My, the moment corresponding to initial yield in the absence of axial thrust, 
and the curvature value by dy, the curvature corresponding to My- 

To have these results in a form suitable for direct use in Eq. 5 it is nec- 
essary to further develop them by obtaining tangent values at various points 
along each curve. These tangent values will be of the form 


* Generalizing, it can be shown that if the member is pin-ended, regardless of 
lateral loads, inelastic or elastic action, etc., 


= P (Ay) . 


For the case of restrained members, a factor taking into account the end 
condition will need to be added to the right hand side of the equation. 
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which can be considered as a measure of the desired internal stiffness term, 
AM/4A®@ |int, in Eq. 5. The tangent curves obtained from Figs. 11 and 12 are 
shown in Figs. 13 and 14. In each case the tangent values are plotted as a 
function of the two coordinates M/My and ¢/¢y. 

If a cosine deflection curve is assumed and the external rate of change 
substituted in 


In these equations, P, is the Euler buckling load for the member of length L. 
APPLICATIONS 


For any member subjected to both end thrust and bending moments, the 
total moment at any section is composed of two parts: one is due to the axial 
thrust times the deflection; the other is independent of deflection (e.g. end 
moments, lateral loads, etc.). That part of the centerline moment which is 
independent of deflection is noted throughout the remainder of this discussion 
as "Mo." 

If it is assumed that the moment producing loads are symmetrically 
placed on the structure and that the deflection curve can be approximated by 
the equation 


= Yo COS (11) 


the external rate of increase in centerline moment with respect to curvature 
will be given by the equation 


4M| 

ext 
However, to have this in the non-dimensional form suitable for use in Figs. 
13 and 14, it is necessary to multiply both sides of the equation by 1/EIL. 
This then, as was shown in the last section, gives 


4M oy 
My 
A M M 
int 
then, 
“9 
° My int Pe 
since 
4M ¢y | _ PL? dy PL? P 
2 2 
ext” * EI Pe 
(12) 
MY lint * 
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where P, is the Euler buckling load for the pin-ended member. 
Since a cosine deflection curve (Eq. 11) was assumed, the centerline de- 
flection can be expressed as 


(Ho 


Substituting Eq. (11a) in the following expression for m, 


(13) 
5» 
= =. (P_) (14) 
My My Pe 


It should here be noted in the derivation of Eq. 14 it was not necessary to 
consider the explicit form of the moment producing loads; only the axial 
thrust, an assumed deflection curve and the length of the member was used. 
Eq. (14) then is the solution to the general inelastic instability problem for 
members, which due to applied loads, deform symmetrically according to 
Eq. (11). 

To illustrate the procedure, consider the following numerical example: 

Given an 8WF31 section, L/r = 80, subjected to a constant axial 

thrust value of P/P, = 0.4; determine the critical values of various 

"moment producing, loading components." (See Fig. 15,). 

For this problem, according to Eq. (12), 4M/A@- A, /My lint = (0.346), 
which from Figs. 13 and 14 give the critical centerline moment* and corre- 
sponding curvature” values as 


Weak Axis 
Bending 
(Rectangle) 


Strong Axis 
Bending 


Substitution of these quantities in Eq. 14 gives the following values for 


m,/My 


*These values have been circled in Figs. 13 and 14. 
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m M 
: 
| 
\ 
| 4m 
AO 
My 
| Mo 
— 0.624 0.904 
My 
0.346 
| 
0.653 1.105 
y 


Strong Axis Weak Axis 


Bending 
Bending Rectangle) 


Since in Fig. 15 the moment at the center of the column exclusive of thrust 
times deflection is equal to 


the critical value of the moment producing loads will be given by the follow- 
ing expressions: 
Me + RL/4 
(2) = (Me + RL/4) cr = 0.398 (for strong axis bending) 
My My 
sal = 0.522 (for weak axis bending) 


Had other loading conditions been imposed, a similar procedure would have 

been used. The expressions for a few typical problems are shown in Table IL. 
In each of these problems, the quantity of primary importance is Mo/My. 

To facilitate solution of other problems where length and P/ Py values are 


different from this one, values of m,/ wy have been obtained for both strong 


axis and weak axis bending in a range of slenderness from 0 S< L/r £120. The 
plots resulting from these calculations are shown in Figs. 16 and 17. They 
are applicable only to those cases of load which produce single curvature 
deformation that can be approximated by a "full" cosine curve. 
Figs. 16 and 17 can be used to solve three basically different types of 
beam-column problems. 
1. Given a constant axial thrust, determine the maximum lateral load 
and/or end moment the member can carry. 
2. Given the lateral loads and/or end moments, determine the maxi- 
mum axial thrust to which the member can be subjected. 
3. Given an eccentricity, determine the maximum axial thrust. 
The first two are solvable directly from the figures. For the third case, 
however, it is first necessary to "lay-off" a line which described the rela- 
tionship between m, and P.* The point of intersection of this "e-line” with 
the given slenderness ratio will determine the critical axial thrust. 


DISCUSSION AND SUMMARY 


In the foregoing presentation a statical method was used for determining 
conditions of indifference. This method implies that an equilibrium position 
is known and it is desired to determine if that position is stable or unstable. 
Solution was achieved by comparing the rate of change in external moment to 


*This line would correspond to the equation mp = P., or 
(mo/My) = (P/Py) (ec/r2) . 
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My 
0-598 0.522 
y 
cr 
=, * (M, + a 


the rate of change in internal resistance moment for a given virtual displace- 
ment of the loaded specimen. If the internal moment increased at a rate 
greater than the external moment, the specimen was stable. If the rate of 
change in external was the larger, the member was unstable. For the case 
where the rates of change were equal, a condition of indifference (incipient 
instability) existed. 


a. Influence of Assumed Deflected Configuration 
For ease of solution, it was assumed throughout this report that the over- 
all behavior of the member could be predicted with sufficient accuracy by 
assuming a deflection curve, the magnitude of which would be governed by the 
strain distribution pattern at the most highly deformed cross-section. This 
simplification eliminated the need of minimizing integrals and resulted in a 
direct determination of an end moment, axial thrust and critical length cor- 
responding to a given stress distribution pattern. However, when using such 
a procedure, a certain amount of error is introduced, the member essentially 
being forced into an unnatural configuration. To determine the magnitude of 
this error, the following three comparisons will be made. 
a. Comparison of strength results obtained for several different 
curves. 
b. Comparison of the deflection and strength results of each of 
these curves with individual solutions obtained by numerical inte- 
gration, and 

c. Comparison with previous published solutions; namely, 
Chwalla's, (1) and Timoshenko's. (14) 


In the first comparison, the three types of configuration shown in Fig. 18 
are considered. The first of these, that of a "full-cosine” curve, was used in 
the development of the general method of solution (Section II} and also in the 
preceding section on "Applications." The second is a parabola, and the third 
is a “partial cosine” curve of the form suggested by Westergaard and 
Osgood.(16) In the first two cases the moment boundary condition is violated. 
Since each of these assumed curves involves one parameter, yo, the pre- 
viously derived instability condition, 


can be directly applied. For the first two cases, 
a. "Full Cosine” curve 


2 
Ad 
int 
b. "Parabolic" curve 

2 

AM; _pl* (16) 

A¢ 

int 


For the "Partial Cosine" curve (assumption "c"), it is necessary to further 
develop AM/A9| ext So that the ends of the "L” length are held at a given de- 
flection, e, rather than maintaining a hypothetical ")" - length member. (See 
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Fig. 19). If such a procedure is carried out a complex equation containing 


functions of Ay, and Afand cosines thereof is obtained which cannot be solved 
explicitly for Ay, in terms of Ag. However, by considering only the first 
two terms of a cosine series, the expression can be reduced and results in 
the following equation: 


2 


2 


Even though this is the same external rate of moment increase as found for 
case "a," the solutions for members "a" and "c" are different, the Me values 
being determined by different equations. 

a. "Full Cosine” curve 


L2 
Me = My - P $9 —>- (18) 


b. "Partial Cosine” curve 


§ is determined by considering the full value of M, and ¢, at the center 
section, i.e. 


(20) 


Mo/My 


The comparison between these solutions is illustrated in Figs. 20 and 21. 
Here is presented a series of interaction curves for various length ratios for 
the 8WF31 member bent about both its strong and weak axes. The "full 
cosine" deflection curve gave the highest strength values. The "parabola" 
gave the lowest prediction. 

A recent paper by John Clark(4) discusses, in part, this same question. 
Therein he not only investigates the applicability of the "partial cosine" 
curve but also a much more complicated three parameter deflection equation. 
Comparison with test results is also shown in that report. In general, how- 
ever, his results indicate that the increase in accuracy obtained by using the 
more complicated deflection expression does not warrant the added work 
required to obtain a solution. 

The second means of comparison, numerical integration, was used as a 
check in a previous report(9) and one figure based on that work is here in- 
cluded to indicate the observed trends. (See Fig. 22) There is difference be- 
tween the load deflection curves obtained by the various methods; however, 
the ultimate carrying capacities of the members are in comparatively good 
agreement for these examples. 

The third comparison is that between previous published solutions and the 
curves of ultimate strength as predicted by the "partial cosine" deflection 
curves presented earlier in this discussion. Since previous information re- 
garding WF section behavior is lacking, comparison is made only for the 
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rectangular section; furthermore, only for two slenderness ratios. 

Chwalla's solution is based on a numerical integration procedure and 
should therefore be considered as being "exact." His derivation is based on 
(a) a typical German steel with a stress-strain curve having a proportional 
limit of 27 ksi, and a yield stress level of 34 ksi. Jezek's is an approxima- 
tion based on (a) a "full cosine” deflection curve and (b) a stress-strain 
curve of the idealized elastic type (see Fig. 1). Timoshenko's solution also 
considers (a) the "full cosine” curve but (b) the stress-strain curve has a 
proportional limit 30 ksi, a yield level of 36 ksi and an ultimate strength of 
56 ksi. 

Fig. 23 shows the comparison obtained between each of these previous 
solutions and the "partial cosine” solution of this report. All curves have 
been adjusted to a yield stress level of 40 ksi. Length adjustment was ac- 
cording to the equation 


ad 
40.0 

In general, Chwalla's solution gives the lowest strength values; Jezek's gives 
the highest. The partial cosine solution is between these two limits, as is 
Timoshenko's solution. 


b. Elastically Restrained Members 

In actual structures, columns are usually restrained at their ends rather 
than being pin-ended, and a recent paper (2) discusses such a condition. In 
that study a method was developed for reducing the actual restrained system 
to one of an equivalent, pin-ended, eccentrically loaded column. Having de- 
termined by that method the equivalent member for a restrained column 
problem, the solution and curves presented in this paper directly apply. 


c. Residual Stresses 

To include the variable of residual stress, it is necessary to develop a 
complementary set of curves for the desired pattern and magnitude of re- 
sidual stress. This development would start from a given M-¢ curve, (9) 
which included the effect of the residual stress pattern, and proceed along the 
same lines previously described. 


d. Load-Deflection History 

It should be pointed out that no special consideration has been given to the 
manner in which a condition of instability is approached; that is, the load- 
deflection history of the member in question. The presented curves and nu- 
merical work assume that an axial load, smaller than that required for fail- 
ure, is first applied and that bending is then increased from zero to its maxi- 
mum value. Under this condition there will be no unloading of previously 
yielded fibers. The resulting critical loads then will be comparable to the 
reduced modulus load for axially loaded members. 

As was clearly demonstrated, first by Shanley(13) and later by Duberg and 
Wilder(5) and others, if an axially loaded column is considered as an imper- 
fect compression member with the imperfection approaching zero, the 
column will start to bend at the tangent modulus load. Increased thrust be- 
yond this value will be accompanied by an increasing rate of lateral deflec- 
tion, the member reaching its maximum strength at some point prior to at- 
tainment of the reduced modulus load. 
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Since one limit of the generalized beam-column problem is the axially 
loaded column, there will be a certain amount of error involved in the method 
of solution developed in this report. While it is felt that for the assumed 
stress strain curve of Fig. 1, the resulting error is quite small, further ex- 
tension of this work to materials having curved stress-strain relations 
(either as a basic material property or as a result of residual stress, etc.) 
should be accompanied by suitable check calculations to indicate the amount 
of error involved in this type of a procedure. The error will be primarily a 
function of the initial eccentricity. Other important factors would be the type 
of cross-section and the stress-strain relation. 

A recent NACA report (17) provided one source of information for these 
check calculations. Therein load-deflection curves are presented (for an 
idealized WF section column) for various initial eccentricities and stress- 
strain curves. 

Since a majority of the practical beam-column problems encountered in 
the analysis of Civil Engineering type structures contain moment producing 
loads greater than that which could be attributed to accidental imperfections, 
this "tangent modulus influence” will more than likely not be a real problem. 

It is interesting to note that the early numerical work on collapse strength 
of rectangular section, beam-columns carried out by Chwalla was based on a 
typical "German Steel” type stress-strain relation (linear up to the propor- 
tional limit, parabolic from the proportional limit to the yield stress level 
and then horizontal to the point at which strain hardening commences). Re- 
cently (2) this work has been corrected and it was found that adjustment was 
needed only when the ratio of the centerline deflection at maximum load to 
the depth of the section was less than 0.1, a comparatively small deflection 
value. 


e. Local and Lateral Instability 
Neither local nor lateral types of instability were considered in the de- 
velopment presented in this paper. The considered condition of instability in 

each case has been that of excessive bending in the plane of the applied mo- 
ment. However, when bending is forced about the major axis of a section 
(that condition most frequently encountered in practice); when there is ap- 
preciable difference in bending stiffness between the two principal axes of the 
section; and when there is little if any lateral support along the member, the 
mode of failure will more than likely be that due to a combined bending and 
twisting action.(8) In practice, however, this will not necessarily be the case 
since members will be frequently restrained by bracing, walls, floors, etc. 
Where sufficient lateral restraint is provided, the problem of bending plus 
twist is nonexistent, and collapse can be determined by the methods previ- 
ously described. 

Another entirely different type of problem can develop when the cross- 
section is composed of connected, relatively thin, projecting elements. For 
those types of members, the projecting elements may individually or collec- 
tively deform in certain local regions and thereby change the original cross- 
sectional properties of the member. 

Assuming the validity of any of the currently available inelastic local buck- 
ling theories, the maximum Strain, € max, that a flange can sustain and re- 
main in the straight position can be predicted. From this can be determined 
the critical value of € max/ y: Using the procedure outlined in Ref. 9 and 
noting that 
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(see Fig. 24), the corresponding critical value of %, /Py for available, rolled, 
wide-flange shapes can be determined and is found to - comparatively large 
(¢o/@y > 3.0). From Fig. 13 it can be seen that such a condition will be 
realized only at extremely low values of thrust, (P,,/Pe < 0.05). 


SUMMARY 


This report presents a virtual displacement solution to the inelastic in- 
stability problem, the method being equally adaptable to elastic, inelastic, 
concentrically loaded, eccentrically loaded and or laterally loaded columns. 

The essence of the solution is that of applying a virtual curvature change 
to the loaded member and equating the rate of change in external moment at 
the most highly stressed section to the internal stiffness resulting from that 
curvature change. For ease of solution, deflection curves of one parameter 


were assumed causing the system to reduce to collocation at one point. It was 


shown that the ultimate strength of such members is relatively unaffected by 
choice of deflection curve, the "full cosine" giving the highest strength pre- 
diction and the parabola giving the lowest prediction. (Figs. 20 and 21) The 
"partial cosine” solution compared favorably with previously available solu- 
tions. (Fig. 23) 

To facilitate solution of problems other than those discussed herein, 


curves have been presented which ensure rapid determination of critical load- 


ing values (Figs. 20 and 21) for any symmetrically deformed pin-ended 
member, 

The presented solutions of this paper apply only to mild structural steel 
and in the strickest sense only to the case where an axial thrust smaller than 


that required for failure is first applied and then held constant as bending mo- 


ments are increased to collapse. This condition, however, can readily be ap- 
proached in Civil Engineering Type Structures. 


ACKNOWLEDGMENTS 


The author wishes to express his sincere appreciation to Bruno 
Thurlimann, who offered many valuable suggestions throughout the develop- 
ment of this approach to the inelastic instability problem. Grateful acknowl- 
edgement is also extended to Paul C. Paris for criticism and suggestions of 
the manuscript. This work has been carried out as part of the project 
"Welded Continuous Frames and Their Components” being conducted under 
the general direction of Lynn S, Beedle. It is sponsored jointly by the Weld- 
ing Research Council and the Navy Department with funds furnished by the 
following: American Institute of Steel Construction, American Iron and Steel 
Institute, Column Research Council (Advisory), Office of Naval Research 
(Contract No. 39303), Bureau of Ships and Bureau of Yards and Docks. The 
helpful criticisms of members of the Welding Research Council Lehigh Proj- 


ect Sub-committee (T. R. Higgins, chairman), and the Column Research Coun- 


cil Research Committee D (N. M. Newmark, chairman and N. J. Hoff, co- 
chairman) are sincerely appreciated. This program is being carried out in 


the Fritz Engineering Laboratory, Lehigh University of which William J. Eney 


is Director. 


692-14 


a 
= 


NOMENCLATURE 


Area of cross-section 


Width of cross-section 


Depth of cross-section 


End eccentricity of axial load application 


Young's modulus of elasticity 


Moment of inertia of total cross-section about the principle axis 
of section under consideration 


Length of member (distance between pins) 


A Length of hypothetical member (used in "Partial Cosine” deflec- 
tion curve assumption) 1 
Ak Change in length of hypothetical member due to a virtual change 


in centerline curvature 


Bending moment 


Mint Internal bending moment due to stress distribution across section 


Mext External bending moment due to lateral loads, end moment, ini- 
tial curvature, and axial thrust times deflection 


Bending moment at which yield point is reached in flexure 


M, Bending moment at the centerline section of the member 
Mo 


That part of the centerline moment which is independent of 
deflection 


End applied bending moment 


4M, Change in value of centerline moment due to a virtual change in 
curvature at that same section 


Axial thrust 


Axial thrust corresponding to yield stress level over entire 
section 


Euler buckling load 
Por Maximum thrust a member can carry 
Concentrated lateral load 
Flange thickness (wide-flange section) 


u Deflection parallel to web in undeformed position ; 
v Deflection perpendicular to web in undeformed position 
y Deflection in plane of moment application 


Centerline deflection 


Change in centerline deflection due to a virtual change in center- 
line curvature 


Ratio of the depth of yield penetration in compression to total 
depth of section 
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Unit strain 


ty Unit strain corresponding to initial yield 
€max Maximum average unit strain on compression flange of wide- 
flange section (strong axis bending) 
Unit normal stress 
Dy Lower yield point stress 
Oy Variable normal stress on tension side of cross-section (positive 
when in tension) 
Curvature 
%o Centerline curvature 
Ago Virtual change in centerline curvature 
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FIG.!| IDEALIZED STRESS-STRAIN 
RELATION 


FIG.2 TYPICAL LOAD- 

DEFLECTION RELATION 

FOR AXIALLY LOADED 
COLUMNS 
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ypical inelastic Behavior 


FIG 3 LOAD-DEFLECTION RELATION FOR ECCENTRICALLY 
LOADED COLUMNS” 
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Cross-Section Stress Distribution 
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FIG.4 BASIC INTERNAL STRESS- 


LOAD RELATIONS FIG TYPICAL BEAM-COLUMN 
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FIG6 CHANGE IN STRESS 
DISTRIBUTION AND RESISTING 
MOMENT DUE TO A VIRTUAL 
CHANGE IN CURVATURE 
(ELASTIC-AXIAL LOAD) 


FIG 7 TYPICAL BENDING RESISTANCE 
RELATION FOR AXIALLY LOADED, 
ELASTIC COLUMNS 
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FIG 8 TYPICAL BEAM-COLUMN 
SUBJECTED TO VIRTUAL 
DISPLACEMENT 


FIG.9 CHANGE IN STRESS DISTRIBUTION 

AND RESISTING MOMENT DUE TO A 

VIRTUAL CHANGE IN CURVATURE 
(ELASTIC- ECCENTRIC LOAD) 


Typical inelastic Behavior 
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FIG.10 TYPICAL BENDING RESISTANCE RELATION FOR 
ECCENTRICALLY LOADED COLUMNS 
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FIG. I! MOMENT-CURVATURE RELATIONS FOR THE 
8 WF SECTION, STRONG AXIS BENDING 
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FIG. 12 MOMENT-CURVATURE RELATIONS FOR 
THE RECTANGULAR SECTION’ 
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INTERNAL RATE OF CHANGE OF MOMENT VERSUS 
MOMENT AND CURVATURE FOR CONSTANT AXIAL 
THRUST VALUES 
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FIG. BEAM-COLUMN 
(used in example problem ) 
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FIG. 17 INTERACTION CURVES FOR 8 WF 3! SECTION 
(WEAK AXIS BENDING) 
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FIG. 20 


FIGI8 ASSUMED DEFLECTION CURVES 
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FIG. 9 VIRTUAL DEFLECTION FOR PARTIAL 
COSINE” ASSUMPTION 
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FIG 23 INTERACTION CURVES SHOWING COMPARISON 
WITH PREVIOUS SOLUTIONS 


FIG.24 STRESS DISTRIBUTION PATTERN WHEN 
STRAINS EXCEED THE ELASTIC LIMIT 
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